Time-Periodic Solutions of the Einstein's Field Equations 
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In this paper, we develop a new method to find the exact solutions of the Einstein's field equations 
by using which we construct time-periodic solutions. The singularities of the time-periodic solutions 
are investigated and some new physical phenomena, such as the time-periodic event horizon, are 
found. The applications of these solutions in modern cosmology and general relativity are expected. 
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1. Introduction. The Einstein's field equations are the 
fundamental equations in general relativity and cosmol- 
ogy. The general version of the gravitational field equa- 
tions or the Einstein's field equations read 
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(1) 



where g^ (fi,v = 0,1,2,3) is the unknown Lorentzian 
metric, is the Ricci curvature tensor, R = g^R^ is 
the scalar curvature, where g^ is the inverse of g^ v , A 
is the cosmological constant, G stands for the Newton's 
gravitational constant, c is the velocity of the light and 
T^ v is the energy- momentum tensor. In a vacuum, i.e., 
in regions of space-time in which = 0, the Einstein's 
field equations (1) reduce to 
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Rp,v — -^g^vR + A<7 M „ — 0, 



or equivalently, 



Rij,u = A.g ( 



(2) 



(3) 



In particular, if the cosmological constant A vanishes, 
i.e., A = 0, then the equation (2) becomes 



R 



1 



-g^R = 0, 



or equivalently, 



Ruv — 0. 



(4) 



(5) 



Each of the equations (2)-(5) can be called the vacuum 
Einstein's field equations. 

The mathematical study on the Einstein's field equa- 
tions includes, roughly speaking, the following two as- 
pects: (i) establishing the well-posedness theory of 
solutions; (ii) finding exact solutions with physical 
background. Up to now, very few results on the 
well-posedness for the Einstein's field equations have 
been established. In their classical monograph 
Christodoulou and Klainerman proved the global non- 
linear stability of the Minkowski space for the vacuum 
Einstein's field equations, i.e. they showed the nonlin- 
ear stability of the trivial solution of the vacuum Ein- 
stein's field equations. Recently, by using wave coordi- 
nates, Lindblad and Rodnianski gave a simpler proof (see 



[jjj]). In her Ph.D. thesis [T^], Zipser generalizes the re- 
sult of Christodoulou and Klainerman [J] to the Einstein- 
Maxwell equations. As for finding exact solutions, many 
works have been done and many interesting results have 
been obtained (see, e.g., 0], [lj|, @)- In what follows, 
we will briefly recall some basic facts about the exact 
solutions of the Einstein's field equations. 

The exact solutions are very helpful to understand the 
theory of general relativity and the universe. The typical 
examples are the Schwarzschild solution and the Kerr 
solution (see Q). These solutions provide two impor- 
tant physical space-times: the Schwarzschild solution de- 
scribes a stationary, spherically symmetric and asymp- 
totically fiat space-time, while the Kerr solution provides 
a stationary, axisymmetric and asymptotically flat space- 
time. 

The study on exact solutions of the Einstein's field 
equations has a long history. In December 1915, Karl 
Schwarzschild discovered the first non-trivial solution to 
the vacuum Einstein's field equations which is a static 
solution with zero angular momentum (see [IB]). The 
unique two-parameter family of solutions which describes 
the space-time around black holes is the Kerr family dis- 
covered by Roy Partrick Kerr in July 1963 (see [101]). 
These solutions are very important in studying black 
holes in the Nature which is just the study of these so- 
lutions (see Q). Various generalizations of the Kerr so- 
lution have been done (see, e.g., [l(| and Q). Gowdy 
@-[3] constructed a new kind of solutions of the vacuum 
Einstein's equations, these solutions provide a new type 
of cosmological model. This model describes a closed in- 
homogeneous universe, space sections of these universes 
have either the three-sphere topology S 3 or the worm- 
hole (hypertorus) topology S 1 ® S 2 . Recently, Ori [l3| 
presented a class of curved-spacetime vacuum solutions 
which develop closed timelike curves at some particular 
moment, and used these vacuum solutions to construct a 
time-machine model. The Ori model is regular, asymp- 
totically flat, and topologically trivial. 

From the above discussions, we see that the exact solu- 
tions play a crucial role in general relativity and cosmol- 
ogy, so it is always interesting to find new exact solutions 
for the Einstein's equations. Although many interesting 
and important solutions have been obtained, there are 
still many fundamental but open problems. One inter- 
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esting open problem is if there exists a "time-periodic" 
solution to the Einstein's field equations. One of the main 
results in this paper is a solution of this problem. 

In this paper, we focus on finding the exact solutions of 
the vacuum Einstein's field equations (3) and (5). We will 
present a new method to find exact solutions. Using this 
method we can construct interesting and important exact 
solutions, for example, the time-periodic solution of the 
vacuum Einstein's field equations. We analyze the sin- 
gularities of time-periodic solutions and investigate some 
new physical phenomena enjoyed by these new space- 
times. We find that the new time-periodic solutions have 
time-periodic event horizon, which is a new phenomenon 
in the space-time geometry. The applications of these 
solutions and their new properties in modern cosmology 
and general relativity may be expected. 

More precisely, our time-periodic solution to the vac- 
uum Einstein's field equations in the spherical coordi- 
nates (t, r, 0, ip) can be written in the following form 



ds 2 = (dt,dr,dO,d<p)(rj llu )(dt,dr,d0,dtp) J 



(6) 



where 



in which 
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' G= 1 + 2stt + sm9 cos (t-r), 
K = r + m In \r — m\ + e sin (t — r), 
M = 0+ sin 0, (8) 

Q = -i(l + 2siu0)fi~. 



In the above, e S (— 4, i) and m€l are two parameters, 



and f2 ± are defined by 

fi ± = |tan0/2|3 ± |tan0/2|- 



(9) 



the Einstein's field equations, one key point is to choose 
a suitable coordinate system. A good coordinate system 
can simplify the equations and make them easier to solve. 
In the study on the Einstein's field equations, there are 
three famous coordinate systems: harmonic coordinates, 
wave coordinates and the Gaussian coordinates. In this 
paper, we consider the metric of the following form 



-K 2 sin 2 9 J 
(7) 



In Section 5, we analyze the singularity behaviors and 
find some new physical phenomena. 

According to the authors' knowledge, (6) gives the 
first time-periodic solution to the Einstein's field equa- 
tions. Here we would like to point out that, by using our 
method, we can re-derive almost all known exact solu- 
tions to the vacuum Einstein's field equations, for exam- 
ples, Godel's solution Jf}, Khan-Penrose's solution [ill ]. 
Gowdy's solution etc. Our method can also be 

used to find exact solutions of the Einstein's field equa- 
tions in higher dimensions which will be of interests in 
string theory. 

2. Lorentzian metrics. The Einstein's field equations 
are a second order global hyperbolic system of highly 
nonlinear partial differential equations with respect to 
the Lorentzian metric g^ (/i, v — 0,1,2,3). To solve 
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where 
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are smooth functions 
) and satisfy go 
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of the coordinates 
In the coordinates 



* tjO * 3j 



'), the line element reads ds 2 = g iiv dx il dx v . 



For simplicity of notations, we denote the coordinates 



, x ) by (t,x,y,z) and rewrite (10) as 
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I u v p q 

v w 

p p 

V q cr 



(11) 



\ where u,v,p,q,w, p and a are smooth functions of the 
coordinates (t, x, y, z). It is easy to verify that the deter- 
minant of {g^v) is given by 
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g = det(<7 M „) = uwp a — v 2 p a — p 2 wa — q 2 wp (12) 
Throughout this paper, we assume that 



g<0 



(H) 



On the other hand, it is easy to see that at least two of 
the functions w, p, a have the same sign. Without loss of 
generality, we may suppose that p and a keep the same 
sign, for example, 



p < and a < 0. 



(13) 



We can easily prove the following theorem. 
Theorem 1 Under the assumptions (H) and (13), the 
metric (g^) is Lorentzian. 

We now consider the solutions of the Einstein's field 
equations. By Bianchi identities, we believe that the gen- 
eral form of the solutions of the Einstein's field equations 
takes one of the following forms 
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(Type I) 
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or 



/ u v p 

w w 

p p 

V a . 



(Type I) 



For type I, the assumption (H) is equivalent to v ^ 0. 
Therefore, by Theorem 1 we have 
Conclusion 1 If 

p<0, a < and v ^ 0, (14) 

tften i/ie metric (r]^ v ) is Lorentzian. 

For type I, we have 
Conclusion 2 If w, p, a are all negative functions and 
v 2 + p 2 + q 2 ^ 0, then the hypotheses (H) and (13) are 
satisfied, and the metric (r]^ v ) is Lorentzian. 

Similarly, for type M we have 
Conclusion 3 If u is positive and w,p,a are negative, 
then the hypotheses (H) and (13) are satisfied, and the 
metric (r]^) is Lorentzian. 

We are interested in finding exact solutions of the Ein- 
stein's field equations of the above types I-IH. 

3. New method to find exact solutions. In order 
to illustrate our method, as an example, we use the 
Lorentzian metric type I to construct some interesting 
exact solutions, in particular time-periodic solution for 
the vacuum Einstein's field equations (2). 

For type I metrics, our method can be described by 
the following algorithm: 




G 2 = -Ap 



G 03 = -Aq 



Goo = — Au 



FIG. 1: The algorithm to construct the exact solutions 



In fact, the equations (2) can be rewritten as 



(15) 



where is the Einstein tensor. Noting the special form 
of type I, we have 
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One of the equations (15) reads 

Gn = -Ar?n. 
Solving the ODE (17) gives 



(16) 
(17) 



I f Pxx . °~xx 1 / Px \ 1 / &X \ 2 P°~ 

, = , exp|y \- 7 + — -- 2 {j) -2( 7 )j R 

(18) 

where v = v a (t,y,z) is an integral function depending 
on t, y and z, provided that (pa) x =^ 0. 
In particular, by taking the ansatz 

p = p(t, y, z) exp{2/(i, x)}, a = a(t, y, z) exp{2f(t, x)}, 

(19) 

(18) becomes 



dx 



v = v f x e f . 



(20) 



According to the algorithm shown in Fig 1, in a similar 
way we can solve other equations in (2), and then we can 
construct many exact solutions of (2). 
Remark 1 Our method can be used to solve the Ein- 
stein's field equations with physically relevant energy- 
momentum tensors, e.g., the tensor for perfect fluid: 
Tjs = (p + p)u 7 us + pg-fS, where p > is the density, 
p is the pressure, u stands for the space-time velocity of 
the fluid with u 7 u 7 = —1. 

Remark 2 The method presented in this paper can also 
be used to solve the Einstein's field equations in higher 
space-time dimensions. 

4- Time-periodic solutions. By the method presented 
in the last section, we can construct many new exact 
solutions to the vacuum Einstein's field equations 



Gfxv - 0. 



(21) 



For example, in the coordinates (r,r,6,tp), taking the 
ansatz in (19) as follows 

p = — 1, a = — sin 2 6, f = ln[r + r + esinr], 
one can obtain an interesting solution of the form 



( Voo Voi V02 

r?oi _0 

V02 T]22 
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(22) 
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where 



Voi 

V02 
V22 



1 + 2sQ + sin 9 cos t H 0, 



Proof. Noting e e (— |, g), we have 
r/oo = G = l + 2£f2+sin0cos(t-r) 



0+ sin6> 



m 



n, 



to 



i(l + 2sin0)Q _ [r + r + esinr], 
— [r + r + e sin rl , 



(23) 



= l + 4e 



COS ; 



cos ■ 



f?33 = — [r + t + e sin r] sin 0, 



in which e is a parameter, to ^ is a constant, and 



.99, 
sm - cos - cos(t — r) 



= exp ( T + T ) +m, 



TO 



|tan0/2|2±|tanfl/2|-2. 

(24) 

The solution (22) belongs to the class of type I, also be- 
longs to the class of type IH. 
Making the transformation 



> i-4N^ 

|cos(t — r)\ 
: l-8|e| > 1 



• 3 

sin - 
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0. 



' t = t + r, 
r = to + exp 
= 0, 



On the other hand, by calculations we have 
Mr 



(30) 



t + r 



(25) 



-G + 



G -G + 
Mr 



G- 



r — to 
2Mr 



Then the solution to (21) becomes, in the coordinates 

(t, r,6,ip), 



r — m r — to 

for r 7^ 0, to and 9 ^ 0, n, 

Mr 



(Mr) 2 
(r — to) 1 



<0 (31) 



ds 2 = (dt, dr, d9, d<p)(r)p V )(dt, dr, d9, dip) 1 



(26) 



where 
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' G= 1 + 2en+ sin (9 cos (*-r), 
if = r + to In \r — m\ + e sin (f — r) , 
M = Q+sin6>, 

Q= -hl + 2sm9)fl-. 



In (28), fi ± are defined by 

± = |tan6»/2|5 ± |tan0/2| _ 3. 






-if 2 sin 2 6> 
(27) 



(28) 



and 



(32) 
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G -G+ QK 
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-G + J^ G-^L -Qif 
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Qif 
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-Qif 




-if 2 
-if 2 sin 2 9 



-ifW*^ 2 



(r — to) 2 
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(29) 



An important property of the space-time described by 
(26) is given by the following theorem. 
Theorem 2 When e takes its value in the interval 
(— g, g), i.e., e 6 (— g,g), the solution (26) to the vac- 
uum Einstein's field equations is time-periodic. 



(33) 

for r 7^ 0, to, 7^ 0, 7r and if ^ 0. In next section, 
we will show that 9 — 0, n, r — 0, to and if = are 
not essential singularities, in other words, they are not 
physical singularities. 

The above discussion implies that the variable t is a 
time coordinate. Therefore, it follows from (27) and (28) 
that the Lorentzian metric (26) is indeed a time-periodic 
solution of the vacuum Einstein's field equations. This 
proves Theorem 2. ■ 

Remark 3 The time-periodic solution (26) was first ob- 
tained by us in early 2007. The authors have presented 
this solution in several conferences for the past year. As 
mentioned above, according to the authors' knowledge, 
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this is the first time-periodic solution to the Einstein's 
field equations. 

Remark 4 In July 2006, the first author discussed with 
Dr Y.-Q. Gu about some ideas presented in this paper. 
We noted that for the case of type I metric, the author 
°f 0/ constructed some exact solutions, none of which 
is time-periodic, because the variable t in these solutions 
can not be taken as the time coordinate. 

Direct computations give us the following property of 
the time-periodic solutions. 

Property 1 In the geometry of the space-time (26), it 
holds that 



Combing (28) and (41) gives 

g = — (£l + ) 2 sin 9 [r + mln \r — m\ + esin (t — r 



Noting (24), we obtain from (42) that 



A 



i4 r' 



(r — to) 2 
(42) 



S r =o,m = {(t,r,9,tp)\ r = 0, m}, 



A 



Sk=o = {(t, r, 9, tp)\ r + to In \r — m\ + e sin (i — r) = 0}, 



A 



(43) 



~d9~ 



dK _ 
~dt ~ 

Or 
2 sin^' 



G - 1 d 2 K 



2M 



dt 2 



1 dG 

2M dt ' 



(34) 



5-^-2, -Xn = Q, 35 
2 sm 9 d9 



So=Q,« = {(t,r,8,<p)\0 = 0, tt} 

are singularities for the solution metric (26). 
By a direct calculation, we have 

R a/3yS = and R a()lS = (a, /3,j, 8 = 0, 1, 2, 3). 
This gives 



dG _ _dG d 2 G 
dr dt ' drdt 



d 2 G d 2 G 



dt 2 ' d9dr 



d 2 G 
d9dt 



and 



dK dK 
dt dr 



d 2 K 
dtdr 



d 2 K d 2 K 



dt 2 



d9dr 



, (36) 



d 2 K 

dtd9' 
(37) 



R ltu = (im,v = 0,1,2,3) 



Rll = R a ^ 5 R a ^ s = 



(45) 
(46) 



dQ 3n+ 

— = Qcot9-— — , 
d9 4 sm 9 



dK _ ldG 
Q ~dT ~ 2~d9' 



d 2 G _ dKdQ 

~dlP ~ dt d9 
(38) 



2Qcos9 _ 3 

~ 4 



Q 2 



1 



(n+) 2 {n+y 



3 O" 1 

2Q cot 9 - - — 



4 sin( 



— sm 



-M. 



(39) 



(40) 



The relations given above will play an important role 
in the future study on the geometry of the time-periodic 
space-time (26). 

There are several important questions which deserve 
further study: (a) what is the topological structure of the 
time-periodic space-time (26)? (b) does the space-time 
(26) have a compact Cauchy surface? (c) an important 
point is to consider the structure of the maximal globally 
hyperbolic part of the space-time (26), the question is 
whether this also exhibits time periodicity. Problems (b) 
and (c) were suggested by Andersson [l|. 

5. Singularity behaviors and physical properties. This 
section is devoted to the analysis of singularities and the 
physical properties of the time-periodic solution (26). 

For the metric (27), by a direct calculation, we have 



A 



g = detfa^) = -M Z K" 



(r — to)' 



■ sm 



(41) 



(45) implies that the Lorentzian metric (26) is indeed a 
solution of the vacuum Einstein's field equation (21), and 

(46) implies that this solution does not have any essential 
singularity. 

According to the definition of event horizon (see e.g., 
Wald [13) j it is easy to show that Se = o^ and Sk=o are 
the event horizons of the space-time (26). If (t, r, 9, tp) are 
the polar coordinates, then the S r= o can be regarded as a 
degenerate event horizon. However, S r=m is a new kind 
of singularity which is neither event horizon nor black 
hole. This differs from that in Schwartzschild space-time, 
in which r = corresponds to the black hole and r = to 
corresponds to the event horizon. In other words, the so- 
lution (26) describes an essentially regular space-time, it 
does not contain any essential singularity like black hole. 
It is an interesting topic to see how to cancel these kinds 
of singularities by making some coordinate transforma- 
tion. Therefore, we have 

Property 2 The Lorentzian metric (26) describes a reg- 
ular space-time, this space-time is Riemannian flat in the 
sense of (44): it does not contain any essential singular- 
ity. However it contains some non-essential singularities 
which correspond to event horizons and some other new 
physical phenomena. 

Property 3 The non-essential singularities of the space- 
time (26) consist of three parts S r =o, m , Sk=o and Se = o t7r . 
5e = o,7r ar £ two steady event horizons; and Sk=o is the 
"time-periodic" event horizon; <S r= o,m o/re two new kinds 
of singularities which are neither event horizons nor black 
holes. In particular, if (t, r, 9, tp) are the polar coordi- 
nates, then S r =o can be regarded as a degenerate event 
horizon. 

According to the authors' knowledge, the degenerate 
event horizon and the time-periodic event horizon are two 
new phenomena in the space-time geometry. 
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We next consider the time-periodic event horizon in 
detail. Without loss of generality, we may assume that 
e and m are positive constants. We have two cases: < 
r < to and r > to. 



In this case, 



mlnm 



<0, 



(58) 



Case I: < r < m. 

Let 

f(r; t) = r + mln \m — r \ + esin (t — r). (47) 
For any fixed t G K, it holds that 

f(r;t) — ► — oo as r — > m (48) 

and 

/(r; t) — ► oo as r — > oo. (49) 
At r = 0, we consider 

to In to + £ sin t = 0, 



i.e., 



sin£ 



mlnrn 



It is obvious that (51) has a solution if and only if 
mlnm 



< 1. 



(50) 



(51) 



(52) 



In what follows, we always assume condition (52). 
Therefore, it follows from (51) that 



/(0;t fc ) = 0, 



where 



tk = 2kir + arcsin 



j to In m 1 
m \—-\ 



(53) 



(54) 



in which k G Z. 

We now divide the discussion into two cases. 
Case 1-1: < to < 1. 
In this case, we have 



to In to 



>0, 



and then 



Therefore, 



mlnm 
arcsin < > > 0. 



. to In to . 

tk = 2fc7r + arcsm ^ f (k = 0, 1, 2 • 



Case 1-2: 1 < to. 



(55) 



(56) 



(57) 



and 



{mlnm ] 
> < 0. 



(59) 



Thus, we shall take 



tk = 2(fc + 1)tt + arcsin < — 



,ln? 



(fc = 0,l,2--.). 

(60) 

In both case 1-1 and case 1-2, by fixing k G {0, 1, 2, • • • } 
and noting (48), we see that there exists a maximum r_ G 
[0,to) such that, for any given r G [0,r_], the equation 
for t 



/(r;t)=0 



(61) 



has solutions. When r = r_, we denote the solution by 



t k . It holds that 



/(r_;tfc)=0. 



(62) 



Summarizing the above discussion, we observe that, for 
case I, the time-periodic event horizons are given in Fig. 
2. 




FIG. 2: Time-periodic event horizons for case I 



Case H: r > to. 

Similar to the discussion of case I, in this case the time- 
periodic event horizons are given in Fig. 3. 




m ro r + r 



FIG. 3: Time-periodic event horizons for case H 



In Fig. 3, ro and r+ arc denned in the following way: 
noting (48) and (49) , we see that there exists a minimum 
ro G (to, oo ) and a maximum r + G (to, oo) such that, for 
any given r G [ro,r+], the equation for t 

/(r;t)=0 (63) 

has solutions. In particular, when r = tq (resp. r = r + ), 
we denote the solution by tk (resp. by t£). That is to 
say, it holds that 

/(r ;t fc )=0 and /(r+;t+) = 0. (64) 

Therefore, we have proved the following property. 
Property 4 The non-essential singularities of the space- 
time (26) consists of three parts r = 0, r = m and 
r + to In \r — m\ + esin (t — r) = 0. r — is a degen- 
erate event horizon, r = m is a steady event horizon, 
and r + mln \ r — m\ + esin (t — r) = are the "time- 
periodic" event horizons. Time-periodic event horizons 
form and disappear in finite times, they propagate time- 
periodically. 

On the other hand, by some elementary matrix trans- 
formations, the metric (rj^ v ) can be reduced to 

f M 2 r 2 1 

= diag |G, - — — — , -K 2 , -K 2 sin 2 6 j . 

(65) 

Noting that (28) gives 

) ~ diag { 1 + 2eO+ sin # cos (t — r ) , 

(n*) 2 ™ 2 * r2 _ r w4 

l + 2en+sin6»cos(t-r)' ' /' 

(66) 

In (66), we have made use of the fact that, when r is 
large enough, it holds that K ~ r because of the second 
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equation in (28). (66) implies that the space-time (26) is 
not homogenous and not asymptotically flat, more pre- 
cisely not asymptotically Minkowski, because the first 
two components in (66) depend strongly on the angle 9. 
Therefore, we have 

Property 5 The space-time (26) is not homogenous and 
not asymptotically flat. 

Remark 5 Property 5 perhaps has some new applica- 
tions in cosmology due to the recent WMAP data, since 
the recent WMAP data show that our Universe exists 
anisotropy (see This inhomogenous property of the 

new space-time (26) may provide a way to give an expla- 
nation of this phenomena. 

Summarizing the above discussion gives the following 
theorem. 

Theorem 3 The vacuum Einstein's field equations have 
a time-periodic solution (26), this solution describes a 
regular space-time, which has vanishing Riemann curva- 
ture tensor but is not homogenous and not asymptoti- 
cally flat. This space-time does not contain any essential 
singularity, but contains some non-essential singularities 
which correspond to steady event horizons, time-periodic 
event horizon and some other new physical phenomena. 

6. Summary and discussion. In this paper we de- 
scribe a new method to find exact solutions to the Ein- 
stein's field equations (1). Using our method, we can 
construct some important exact solutions including the 
time-periodic solutions of the vacuum Einstein's field 
equations. We also analyze the singularities of the time- 
periodic solutions and investigate some new physical phe- 
nomena enjoyed by these new space-times. 

We remark that, by using our method, we can obtain 
almost all known solutions to the Einstein's field equa- 
tions. Our method can also be used to find exact solu- 
tions of the higher dimensional Einstein's field equations, 
which play an important role in string theory. The struc- 
tures of these new space-times, the behaviors of their sin- 
gularities and some new nonlinear phenomena appeared 
in the time-periodic solutions are very interesting and 
important. We expect some applications of these new 
phenomena and the time-periodic solutions in modern 
cosmology and general relativity. 
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